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I. Introduction 



^ , Gravitational lensing is an important and effective window to look for signatures of peculiar 

astrophysical objects such as black holes (BH). This field of activity has lately attracted a lot of 
interest among the physics community. Early works focussed on the lensing phenomenon in the 
weak field (for a review, see [1]), but weak field results can not distinguish between various different 
solutions that are asymptotically flat. What one needs for this purpose is a method of calculation 
in the strong field regime. Progress in this direction have been initiated by Fritelli, Kling and 
Newman [2], and by Virbhadra and Ellis [3]. However, the difl[iculty is that, in the strong field, 
light deflection diverges at the photon sphere. (The conditions for the existence of photon surfaces 
have been rigorously analyzed by Claudel, Virbhadra and Ellis [4]). By an analytic approximation 
method, Bozza et al [5] have shown that the nature of divergence of the deflection angle becomes 
logarithmic as the light rays approach the photon sphere of a Schwarzschild BH. This method has 
been successfully applied also in the Reissner-Nordstrom BH [6] . Virbhadra and Ellis [7] have further 
extended the method of strong field lensing to cover the cases of Weak Naked Singularity (WNS) and 
Strong Naked Singularity (SNS). Bozza [8] has subsequently extended his analytic theory to show 
that the logarithmic divergence near the photon sphere is a generic feature for static, spherically 
symmetric spacetimes. This work is a remarkable step forward in the arena of gravitational lensing. 
Bhadra [9] has applied the procedure to BHs in string theory. The extension of the strong field limit 
to Kerr BH has also been worked out recently [10,11]. All these investigations have indeed thrown 
up a richesse of information about the signatures of BH via lensing mechanism. 
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There is however another exciting possibiHty that has not received enough attention to date: It is 
lensing by stellar size traversable wormholes (WH) which are just as interesting objects as BHs are. 
WHs have "handles" (throats) that connect two asymptotically flat regions of spacetime and many 
interesting effects including light propagation, especially in the Morris-Thorne-Yurtsever (MTY) 
WH spacetime, have been extensively investigated in the literature [12]. WHs require exotic matter 
(that is, matter violating at least some of the known energy conditions) for their construction. The 
idea of this kind of matter has received further justification in the notion of "phantom field" or "dark 
matter" invoked to interpret the observed galactic flat rotation curves or the present acceleration of 
the Universe. Some works on lensing on a cosmological scale involving dark matter do exist [13,14] 
but they have nothing to do with Lorcntzian WHs on a stellar scale. Nonetheless, it might be 
noted that recent works by Onemli [15] show that the gravitational lensing by the dual cusps of the 
caustic rings at cosmological distances may provide the tantalizing opportunity to detect Cold Dark 
Matter (CDM) indirectly, and discriminate between axions and weakly interacting massive particles 
(WIMPs). It is also to be noted that local, static WH solutions threaded by phantom matter have 
also been worked out recently [16]. 

Work in the direction of WH lensing has been initiated by Cramer et al [17] not very long ago 
and recently Safonova et al [18] have investigated the problem of lensing by negative mass WHs. 
A most recent work by Tejeiro and Larrahaga [19] shows that Morris-Thorne type WHs generally 
act like convergent lenses. Unfortunately, work on WH lensing, let alone the strong field analysis, is 
still relatively scarce though observables in WH lensing have the potential to serve a dual purpose: 
They would establish not only the WH itself but also throw light on the existence of classical exotic 
matter. This fact provides the basic motivation for the present theoretical investigation. 

We shall investigate the strong field lensing phenomenon in the WH solutions belonging to the 
Einstein minimally coupled scalar field theory (EMS) as well as the brane world model. (It is to 
be noted that novel effects of the scalar field on gravitational lensing have been analyzed in Ref. [20] 
way back in 1998 in the context of point like naked singularity lens.) Importance of the EMS 
theory need not be repeated here. Suffice it to say that it is the simplest scalar field theory. It 
can be connected to the vacuum Brans-Dicke theory via the so-called Dicke transformations and to 
the vacuum heterotic string theory. The static WH solutions in all these theories have been well 
investigated [21-30]. On the other hand, brane theory is a completely different proposition of great 
interest. The brane world paradigm envisages that only gravity propagates in the 5-D bulk while all 
other fields are confined to the 4-D brane. This idea leads to newer models of local self-gravitating 
objects. It would therefore be interesting to calculate the lensing effect in these models, especially 
in the strong field limit. 

Generically, the brane world BHs are far richer in structure than ordinary BHs as they embody 
a synthesis of wormhole and black hole features. That is why we refer to those objects here as 
WH/BHs. For instance, the effective stress energy tensor could violate some of the energy condi- 
tions, though it need not always be the case. This feature is not unexpected as the stress tensor 
contains imprints of the nonlocal free gravitational field existing in the 5-D bulk which contributes 
negative energy [31]. Several observable effects of the extra dimension on quasar luminosity in the 
rotating models have been recently reported [32]. In the context of spherical symmetry, the extra- 
dimensional bulk contribution essentially implies a correction to the Schwarzschild solution but its 
horizon structure remains undisturbed. The brane theory we have in mind is described by the RS2 
framework, that is, a single brane in a Z2-symmetric 5-D asymptotically anti-de Sitter bulk in which 
only gravity propagates while all other fields are confined to the brane [33]. Strong field lensing 
in one of the brane world BHs [34] have been carried out in Ref. [35] (and the weak field lensing is 
calculated in Ref. [36]). Lensing in another class of brane world BH (see below, Sec.IVA) has been 
investigated by Whisker [37] . Authors in Refs. [35] and [37] have shown that such BHs could produce 
observables that arc significantly different from the Schwarzschild BH. 

In this paper, we shall apply the strong field limit procedure, due to Bozza [8], in the standard 
lensing (distinct from retrolensing) phenomenon by static spherically symmetric WH solutions in 
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the EMS theory and by the WH/BH solutions in the brane theory. (Lansing in the weak field regime 
has been investigated in Refs.[38].) We show that more spectacular differences can appear in the 
observables in the strong field limit. This is our key result. 

The paper is organized as follows. In Sec.IIA, wc outline the procedure of the strong field limit 
including the expressions for observables in Sec.IIB. In Sec.IIIA, we deal with the lensing by a 
massive WH and in Sec.IIIB, with the zero mass WH. In Sec.IVA, we point out that the brane 
world BH, considered recently in Ref. [37] can also be interpreted as a self dual WH harboring a 
globally strong naked singularity. Sec.IVB reveals characteristic features of the strong field lensing 
by a different brane world WH/BH solution. In Sec.V, we point out certain important aspects of 
the considered WH solutions that should be useful in understanding the lensing behavior. Finally, 
Sec. VI summarizes the results. 

II. Strong Field Limit 
A. Deflection angle 

We assume that the asymptotically flat spacetime describing a BH or WH is centered at L which 
serves as the lens. The observer O and the source 5, which is to be lensed, are positioned in the flat 
region on cither side of L, but not necessarily along the same line. This is a plane configuration of 
ordinary lensing, as distinct from retrolensing where both O and S are positioned only on one side 
of L. Let / be the location of the image of S as observed by O and that the extended IS segment 
meet the extended OL segment at X. Defining the angles as {{OL,OS) = (3, {{OL,OI) = 6, the 
lens equation follows from the plane geometry [3]: 

tan /3 = tan 61 - [tan 61 + tan (a - 61)] (1) 
Dox 

where Dqx — Dol + Dlx and Dpq is the Euclidean distance between P and Q, and a is the 
deflection angle. 

The generic spherically symmetric static metric for our purposes is (we take SttG = c = 1): 

rfs^ = A{x)dt^ - B{x)dx^ - C{x)dn^ (2) 

where dfl"^ = dO'^ + sin^ Odip^ is the metric on a unit sphere. The photon sphere a solution 

of the equation 

C{x) A{x) ^ ' 

in which the primes represent derivatives with respect to x. The impact parameter u is defined in 
terms of the closest approach distance x = xq as 



The minimum impact parameter is given by 



A{xo) 



IC{xp 



(5) 



y A(^Xps) 

Prom the equation of photon trajectory, it is easy to derive the deflection angle 

a{xo) = -TT + Iixo) (6) 



I{xo) = / , , (7) 
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Bozza's procedure [8] for the strong field limit is based on the following conditions: (a) The 
photon sphere x = Xps must exist, (b) The functions A,B,C,A',C' must be positive for x > Xps, 
(c) There should exist a static limit, or horizon where A{xs) = 0. The last condition is sufficient but 
not necessary. Then, define y = A{x) , yo = ^(a;o) and 



1 -yo 

and rewrite the integral I{xo) as 



^ (8) 



/(xo) = / R{z,xo)fiz,xo)dz (9) 



R{z,xo) = ^^{l-yo)Vc~o (10) 

/(^, xo) = I ^ == (11) 
\Jyo-[0--yo)z + yo\%- 

where all functions without the subscript are evaluated at a; = ^~^[(1 — yo)z + yo]- The function 
R{z, xq) is regular for all values of its arguments, but the function f{z, xo) diverges as 2: ^ and it 
expands to second order like 

f{z, Xo) ~ fo{z, Xo) = ^ , r, n (12) 

where the parameters a\ , /3i depend on the closest approach xq as 

"i = ^K2/o-Co^(,] (13) 
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A = [2CoC^A[,2 + (CoC^ - 2C'^) yoA^, - CoC^yo^^,'] (14) 

Then the integral /(xq) is resolved into a regular and a divergent part and the latter gives the 
deflection angle to order 0{xq — Xps) as 

a(e) = -aln( -1^+6 (15) 

V Ups J 

where 

a=*^ (16) 
28 

b= -TT + bR + aln-^ (17) 

Vps 

bR= g{z, Xps)dz + 0{xo - Xps) (18) 
Jo 

g{z, Xps) = R{z, Xps)f{z, Xps) - R{0, Xps)fo{z, Xps) (19) 

Pps = Pi \xo=Xp, ^Vps = Ai^Xps) (20) 

The function g{z, Xps) is regular at z = [8]. 
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B. Observables 



The relativistic images of the source are greatly demagnified in comparison to weak field images 
because the photon trajectories wind several times around the photon sphere before emerging outside. 
Yet, best results are obtained when the source S, lens L and the observer O are highly aligned. In 
this case wc can assume that the angles 9 and /3 are small, but a = 2n7r + Aa„ , n G Z where Aa„ 
is the residual angle after the trajectories wind the photon sphere n times. Under these conditions, 
the lens Eq.(l) reduces to 

e = (}+^Aan (21) 
Defining a{6^) = 2mr, and using Eq.(15), we can write 

^°=^(l + e„) (22) 

where 

e„ = e^^-^"'^)/" (23) 
The position 9^ and the magnification of the nth relativistic image are: 

il3-9l)Dox 

9n = 9„-\ — (24) 

O-J^LX^OL 



I I e.nul^{l + en)Do x 



Now we bunch all the images together at 9^0 = Ups/DoL, so that the outermost single image 
appears at 9\. Then define the observables 

s = 9i-9oo (26) 
" - (27) 



Eoo 
n=2 

which, respectively, are the separation and flux ratio between the bunch and the outermost image. 
Using the relevant expressions, they simplify to 

s = 9^e^-'^''^l'^ (28) 

r = e^'^/" (29) 

We shall calculate the strong field cocfficicmts a, h and the observables s, r for some physically 
interesting WH solutions in the EMS and the branc world model. 

III. EMS theory 

The field equations of the EMS theory are 

R^, = (30) 



= (31) 

where $ is the minimally coupled scalar field and k is a constant free parameter. Note that the 
above equations are just the conformally rescaled vacuum Brans-Dicke equations [24,30]. Clearly, 
all the results in the sequel can be easily transcribed into those of Brans-Dicke theory and further 
on, to string theory [29]. 
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A. Massive WHs 



A well known class of solutions of the EMS theory is the Janis-Newman-Winnicour ( JNW) [39] 
solution (or a variant of the Wyman [40] solution): 



A{x) 



2™ 

X 



,B{x) 



2m y 
X ) 



,C{x) 



^ _ 2mV ^ 
X ) 



(32) 



$(x) = 



where M is the ADM mass given by 





In 


1- 


2m" 


V 2k 


X 


7 = 


M 
m 







Q 

X 



(33) 
(34) 



= m^ - -nq^ 



(35) 



m 



2(1-7^) 



(36) 



is the asymptotic scalar JNW charge. In the field equations (30), we have introduced a new constant 
parameter k that does not appear in the observables but facilitates the analysis of the nature of the 
EMS solutions. With a positive sign on the right hand side of Eq.(30), the stress tensor represents 
ordinary scalar matter with positive energy density. The solution (32,33) then has a globally strong 
naked singularity at x = 2m when 7 < 1. However, with a negative sign on the right hand side, the 
stress tensor represents energy condition violating exotic matter necessary for constructing WHs. 
Now, this negative sign can be achieved in two ways: (i) Take k = —2 (that is, break all the energy 
conditions "by hand" or assume that this sign comes as an input from another theory) and keep 
$ real or (ii) Take k = 2 but make $ imaginary or which the same thing, q imaginary. The latter 
case also throws up a negative sign on the right side of Eq.(30) and is completely physically valid 
as discussed by Armendariz-Picon [41]. In either case, he solution represents the spacetime of a 
symmetric traversable wormholc [22,24,27,28]. That there are two asymptotic regions can be best 
seen by transforming the metric (30) into isotropic coordinates via a radial transformation 



x = pll 



m 



(37) 



in which case the solution reduces to the Buchdahl solution [42] of 1959 given by 



A{p) = 



1 rn 

2p 



27 



1 + 



2p. 



• B{p) 



^ m 

~Yp 



2(1-7) 



1 + ,C{P)=P'B{P) (38) 



2(1-7^) 



In 



2p 



1 + 



2p 



£ 
P 



(39) 



The solution is invariant in form under radial coordinate transformation P = and hence one 
asymptotic region occurs at p = 00 and the other at p' = 0, the two coordinate patches meeting at 
p = p' = ^. The WH throat occurs at pth = t (7 + -\/7^ — 1^ and the requirement that pth be real 
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and positive demands that 7 > 1. This is the WH condition. The energy density pn and the scalar 
curvature R for the solutions (38) and (39) become 



PD 



l\m2(l-72)/ m\-2(i+7)/ rn\-2(i-7) 



(^-2) (^0) 



^ „ 24/ m\-2(2+7)/ TO\-2(2-7) 

ii = 2mV(p+2-j (P-^) (41) 

Clearly, /Od < for 7 > 1 so that the Weak Energy Condition (WEC) is violated. For M ^ 0, and 
for the case (i), we have 7 = — j ^'^d fo'^ the case (ii), defining q — iq'with q' > 0, we have 

7 = —^^==. Thus, 7 increases beyond unity if q (or g') is non-zero. The photon sphere appears 
at 

TTi r J n 

Pp, = -[27±V47'-lJ (42) 

It is clear that > pt/i so long as 7 > 1. All the functions A{p),B{p),C{p),A'{p) and C'{p) are 
positive for p > ppg. There is also the so-called static limit at Ps = ^ at which A{ps) = 0. But the 
surface Ps = ^ is a, strong naked singularity. However, pps > pth > Ps for 7 > 1 which implies that 
the photon sphere hides the throat and the naked singularity. The situation resembles the lensing 
scenario by Weakly Naked singularity (WNS) defined by Virbhadra and Ellis [7] to the extent that 
the naked singularity is hidden under the photon sphere. The occurrence of a throat hiding further 
the naked singularity is a new feature in the present case. However, the main difference is that the 
Virbhadra-EUis choice of 7 is still less than unity for g' > M since they defined 7 = , 

The calculation of the strong field limit coefficients becomes awkward in the isotropic coordinates 
and it is more convenient to use the metric (32) which is in standard coordinates. Then the photon 
sphere appears at Xps = m{2'y + 1). Without involving any loss of rigor, all that we need to do is to 
take the WH range of 7 from from the foregoing analysis. Now. in the case of Schwarzschild lensing, 
the value u — Ups = 0.003 involves an error of only 0.4% from the exact position of the outer image 
[8]. Taking this value as the starting point and using u = ODql, the coefficients become 

a = 1 (43) 

h - ^^h , , [(27+1)^ -(27 -1)1^(27+1) 

b - -n + bn + ln 272(27-1)27-1 (^4) 

bR = 0.9496 - 0.1199(7 - 1) + 0(7 - 1)^ (45) 
Up. = S^l±2riL (46) 

2(27-1)^-5 ^ ' 

^ [(27+l)^-(27-l)T .47. 
Pps 472(472-1)7-1 ^^'> 

«(7) = -aln(^]+b (48) 

V Ups J 

It was shown in Rcf. [8] that the deflection angle a{^/) decreases from the Schwarzschild value in 
the range of naked singularity (7 < 1). In contrast, the deflection angle a(7) actually increases 
from the Schwarzschild value with the value of increasing 7 in the WH range (7 > 1) as will be seen 
in the Table I below. This behavior is markedly different from the case of naked singularity or the 
Schwarzschild BH. 

B. Massless WHs 
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Consider a WH for which the ADM mass M = 0. The cases (i) and (ii) mentioned in Sec.IIIA 
respectively give = + and = + g'^ which imply that both m = and q = q' = 0. This 
is a trivial case. But we can also have M = rwy = by putting j = 0, m 0. However, we must 
remember that, physically, the solar system tests fix 7 ^ 1 while one is free to choose to = to 
achieve M = 0. This notwithstanding, we consider the reverse case (7 = 0, m 7^ 0) here only as a 
mathematically interesting possibility. Thus, taking K — 2,we have from Eq.(35) that 

q"^ = (49) 

implying that the gravitational stresses due to to and non-gravitational stresses due to q exactly 
balance each other. This is an extremal situation. Though M = 0, and we should not expect any 
deflection at all, the spacetime is not flat. It is conceptually a classic example of Wheeler's "charge 
without charge" [43], and it is a stable WH [41]. 

Due to our present choice of k, < as argued before, and so rr? < 0. Let us take m — —im'. 
Then, we have q'"^ = to'^ and moreover the energy density pn and scalar curvature R at the throat 
Pth = ^ =\ are given by 

1 2 

^^ = -^;^ = -^ (50) 

It would be interesting to analyze the effect of this massless curvature on the light rays. Also, as the 
radial variable p ^ 00 in Eqs.(40) and (41), both and i? ^ implying that the zero mass WH 
solution is asymptotically flat and also perfectly nonsingular everywhere without a horizon. The 
solutions (38,39) become 

A{p) = 1, B(p) = (^1 + , C{p) = p2 (^1 + , ct,(p) ^ ^ (51) 

and the throat occurs at pth = ^. It also represents the Ellis [44] "drainhole" particle model. The 
Eqs.(51) can be expressed in proper distance Z = p — ^ in a quite familiar form 



= dt^ - dl^ - {f + q'^)d^\ $(/) = ArcTan j (52) 

where we have used the identity ArcTanix) = |ln in $(p) of Eq.(39). The photon sphere 

exists and it appears at Ppg = ^ which coincides with the nonsingular throat radius pth- This is 
an extremal situation. Looking at the necessary conditions, we see that (a) is satisfied even though 
A anywhere. This only points to the fact that (c) is not a necessary prerequisite. The condition 
(b) is marginally satisfied since all the desired functions are positive and nonzero for p > y except 
that A' = 0. [Note incidentally that in the standard coordinates of the metric (1), the throat is 
at Xfh = m{l + 7) for 7 7^ and therefore xth = rn for 7 = 0. The photon sphere occurs also at 
Xps = TO, but the difficulty is that C(x) ^ for x > to in violation of the condition (b)]. Because 
A' = 0, the functions R{z, xq) and i?(0, Xm) diverge, and consequently do the coefficients a and bn- 
Since the object is massless, there is no possibility to consider the Schwarzschild lensing as a starting 
point in the strong field analysis as we did before. 

Perlick [45] has discussed detailed lensing properties of the Ellis drainhole given by the metric 
(52). Ho discussed, in terms of an exact lens equation, the cases that observer and light sources are 
(i) on different sides and (ii) on the same side of the WH's throat. If the observer is closer to the 
throat as the light source, the behavior of the bending angle is similar in both the cases. In terms 
of the metric form (51), we can integrate the deflection angle, Eq.(6), as 

a{p) = -TT + 4ArcTanh (53) 
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where p = po is the closest approach. The function a{p) is real only in the range \p\ < po and this 
simply means that a light ray that starts from one asymptotic end and passes through the throat 
of the WH can not go back to the same asymptotic end from which it has started. It is very well 
possible that the light ray goes to the other asymptotic end, as already discussed by Ellis [44] . [The 
plot of a{p) vs p for po = pth = \ (with units q' = 1) is same as the fig. 8 in Ref.[45] derived 
earlier.] Thus, there are two classes of light rays that start from one asymptotic end: Members of 
the first class turn around before they reach the neck, members of the second class pass through the 
throat and proceed to the other asymptotic end. The borderline cases between the two are light rays 
that asymptotically spiral towards the photon sphere at the throat. These features are qualitatively 
similar to the light trajectories starting at infinity in the Schwarzschild spacetime: Members of 
the first class turn around before they reach the photon sphere, members of the second class pass 
through the photon sphere and proceed to the horizon. The borderline trajectories are those that 
spiral towards the photon sphere. 

It is also possible to calculate the deflection angle, Eq.(6) using the familiar proper form of the 
drainhole metric, Eq.(52), considered in Ref.[45]. The minimum surface area 47rg'^ appears at the 
throat I = which is the same as pps = \. The extremal situation ptt = Pps = ^ now translates 
into Ith = Ips = Taking the closest approach at I = Iq = a (say), we get, for the exact deflection, 
an elliptic function 

2\/l + a^EllipticK [-^ 
a{a) = -TT H L_U. (54) 

where EUipticK (x) is a particular case of hypergeometric function. With units in which q' — I, we 
immediately find that a{a) oo (capture) as a — » 0. The plot a{a) vs a in the range < a < oo 
again shows that the stable massless WH acts qualitatively like a Schwarzschild deflector. What 
is interesting is that the Keplerian mass M is zero, yet light rays coming from the source respond 
to this configuration and images the source. The strong deflection limit around a ~ in Bozza's 
formalism [8] is given by 

a(o) ~ -TT - 2 ln(a - |-) + 2 ln(2g') (55) 

In terms of the distance OL, the same is given by a ~ — tt — 21n(OL) + 21n(4g'). These plots 
approximate the exact deflection pattern perfectly well. 

IV. BH/WHs in the brane theory 

The 5-D Weyl tensor when projected onto the brane produces a trace-free tensor appearing 
in the Shiromizu-Maeda-Sasaki brane field equations [46] 

= -^46; - 4t; - 4k - K (56) 



n;: ^ ^[t;t„^ - tt; - s^ir^^^r^^ - It^)] (57) 
A4 = ^4 (as + 1^^'^ (58) 

kI = SttGn = Ki^/Gn (59) 

where Gat is the Newtonian Gravitational constant (we had earlier put 8t:Gn = 1), A4 and A5 are, 
respectively, the 4-D and 5-D cosmological constants, A is the brane tension. Visser and Wiltshire 
[47] worked out an algorithm for finding solutions when matter fields are present (T^ 7^ 0) on the 
brane. To separate the observable effects of pure bulk gravity from those due to ordinary matter on 
the brane, we set T; = 0. As we are interested in the local self-gravitating objects, we can ignore the 
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cosmological A4 term. The trace of the vacuum brane field equations then simply gives R = where 
R is the Ricci scalar. This equation is solved to derive different classes of brane world BH/WHs. 

The general class of 4-D solutions is given by the metric of the form (2) where A{x) and B{x) are 
two well behaved positive functions for x > Xh and have a simple zero at x = Xh defining the horizon. 
The singularities, if any, of the BH solutions when propagated off the brane into the 5-D bulk may 
make the AdS horizon singular ("black cigar" [48]). However, several classes of nonsingular, static, 
spherically symmetric solutions of the brane world model have been proposed almost simultaneously 
by Germani and Maartens [49] and by Casadio, Fabbri and Mazzacurati (GMCFM) [50] and some 
quantum properties have also been investigated [51]. Under certain assumptions on the behavior of 
the metric functions, Bronnikov, Mclnikov and Dchncn [52] have shown that the generic solutions can 
have Rx S-^ topology of spatial sections. Assuming asymptotic flatness at large x, the global causal 
structure of such solutions coincides with a section of the Kerr-Newman nonextremal solutions. We 
shall specifically consider below two important GMCFM classes (I and II) of solutions. 

A. GMCFM I solution 

In the units such that 2m = 1, the metric components are [49,50]: 

2 _^ 

(60) 

in which k, A are arbitrary constants. Lensing in this spacetime has already been investigated by 
Whisker [37], but some additional observations seem to be in order. This is actually a self-dual 
solution of ii = spacetimes with two asymptotic regions. For different domains of the constants, 
this solution represents a Schwarzschild BH, naked singularity and traversable wormholes [53] . Only 
for X = and A = 1, we have a Schwarzschild BH. In isotropic coordinates x = p{l + j^)"^, the above 
metric becomes 




A{p) = 



1^' 



,i?(p)=(l + ^j ,C{p)=B{p)p' (61) 



The equation A{ps) = gives ps = -But at p = ps, there appears a naked singularity as can 

be seen from the following 

2816xp3 1408xp3 

Pjj=0,Pp = , ,p± = , (62) 

'"^^ (l+4p)6VI(^'^^ {1 + Apf^) ^ ' 

where is the density, Pp, p± are the radial and cross radial components of pressure. The equation 
of state is that of so-called "dark radiation" given by pn — {pp + Ipx) = 0. To get into Whisker's 
notation, one has only to identify A = 1 + e and >f = — e so that ps = ^-g^. Due to the negative sign 
before Pp, Averaged Null Energy Condition (ANEC), which is the weakest, is violated. The metric 
(61) then represents a traversable symmetric WH with the throat occurring at pth = \- In order 
that Pth > Ps, must have e < j. This is the condition for traversability though this condition 
is not strictly needed in the strong field lensing calculation because the light rays are assumed to 
travel only up to the photon sphere and not up to the throat. Since both the throat and the naked 
singularity are hidden below pps, we have here a situation like the EMS case investigated above. 

However, it can be shown [27] that the total amount of ANEC violating matter is Qanec = 
— xln/3 which diverges logarithmically with p. Therefore, unless some technical modifications 
to the solution are made (e.g., as in Ref.[54]), the only way to remove this divergence is to set 
X = which then produces the trivial Schwarzschild BH solution. However, this divergence is not 
a problem as the total gravitating mass is positive and finite as explicit calculations will show in 
Sec.VC. 
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B. GMCFM II solution 



There another solution, described below, that also represents Schwarzschild BH, naked singularity 
and traversable WHs which in the unit 2m = 1 has the form [49,50]. 

A{x) = 1 - \,B{x) = ^-^^^^L_^c{x) = (63) 

The nature of the solution depends on various choices of a constant adjustable parameter o inter- 
preted as a bulk induced "tidal charge" - a Weyl tensor projection from the 5-D bulk into the brane. 
The 4-D effective stress tensor components are: 

4a -3 __4a-3 _ (2a; - l)(4a - 3) 

~ x\^x-2>f'^^ ~ ~a;2(4x-3)'^^ " x\^x-2,f ^ ^ 

For further interesting aspects of this spacctime, sec Scc.V below. The horizon appears at x^ = 1. 
The spacetime structure depends on the parameter r; = a — |. Let us state the various cases [52]: 
(i) If rj < or < a < |, the structure is that of a Schwarzschild BH with a spacelike singularity 
at = |. (ii) If?7>0or| <a<l, then the solution describes a non-singular BH with a WH 
throat at xth = The causal structure is that of the (1 + 1) dimensional subspace of a nonextremal 
Kerr BH solution, (iii) If a = 1, then we have a double horizon at a;/, = 1 with a timelike curvature 
singularity at Xs = | . The global structure is that of an extreme Reissner-Nordstrom BH and finally 
(iv) The range a > 1 corresponds to a symmetric traversable WH with its throat occurring at either 
Xth = 1 or Xth = 0,. For a = |, one recovers the Schwarzschild solution. 
The photon sphere and the minimum impact parameter are given by 



x^s = 77 (65) 



3\/3 



(66) 



which are independent of the tidal charge a. Thus, so long as a < 1 as in the cases (i)-(iii), we 
see that the photon sphere covers the surfaces of event horizon or singularities. In case (iv), too, 
the same situation occurs if 1 < a < |. Thus, all physically meaningful solutions satisfying the 
conditions (a)-(c) of Sec. HA are contained in the range < a < |. We can not take a > | because 

The relevant 



in this case, the WH throat radius Xth = a exceeds that of the photon sphere Xps- 
coefficients work out to: 

3 

ai = 2 



(67) 



/3i 



(68) 



Ax 



ps 



a = 



b = — TT + ha 




((X iCps)(3 Xps) 



Xpsi^ AXps) 



2 y Xps){^ Xps) 

bR= g{z, Xps)dz + O{xo - 
Jo 



In 



6 — 2a;„ 



1 



(69) 

(70) 

(71) 
(72) 
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3V2 



z + l 



zV3 - 2a z Y (3 - 2^) [3 + 2a(z - 1)] 



(73) 



All the quantities above are well defined for a < |. The integral has no divergence on [0, 1] but 
its analytic evaluation in closed form is rather unwieldy. However, we can easily expand g{z,Xps) in 
powers of z: 



, 6i(4a + 3) &i(64a-30a2-33) ^, 



(74) 



(18 -12a) ' 8(3-2a)2 

where 61 = 3\/2(9 — 6a)~i which shows that g{z,Xps) is perfectly regular at 2; = 0. Since the 
solution under consideration resembles that of Schwarzschild in many ways, especially, the photon 
sphere appears exactly at the same value, we can, up to a good accuracy, consider photon orbits for 
u — Ups = 0.003. We can then find the corresponding value of z by employing the expression [8] 



where c is a constant. It turns out that 



c = A 



ps 



1 



C3 2(1 - A)2 

and so xq = 1.554. Prom the definition that z = a,t xq = Xps, we can write 

A{xo) - A{xps) 



z = 



1 A.(^Xps^ 



(75) 



(76) 



(77) 



which gives z = Zmin = 
Schwarzschild value a 

bR 



0.035 corresponding to u — Ups = 0.003. By a Taylor expansion around the 



we now obtain 



Xps) |a=3 dz+ [a 



dg 
da 



dz + Oia 



Therefore 



bR ~ 0.9496 - a 



X 1.565 + a 



(78) 



(79) 



The neglected higher order terms are smaller due to the gradually diminishing factors in the powers 
of (a — I) for < a < |. The deflection a{xo) as a fimction of the closest approach distance xq 
now works out to 



a{xo) = 6i{ - 7r+ ifiln 



(18-6a;o)V3(a;o-l) 



2(a;o-l)v'xg-3(a;o-l) 



where 



n = 



' a:o(3-4a:o) 
(3 - xo){a - xo) 



(80) 



(81) 



Using the Schwarzschild value ^00 = 16.87 jj. arcsec, the expressions for a, b, r , s and u as a 
function of closest approach xo turn out to be 



n 



a = 



b = bR- TT 



In 



6 — 2x0 
Xo-l 



(82) 
(83) 
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s = — 



r = Exp 
33.74(a;o - 3) 



47r 

IT 

Exp 



2{bR - 37r) 



u 



1 



(84) 
(85) 

(86) 



Prom the above expressions, it is evident that, for the tidal charge value a — |, the values for a{xo), 
a,b, r and s differ significantly from other choices of a within the chosen range, especially near the 



photon sphere, xq 



J, the relevant expressions become 



6 -4a 



-77 + 2.123- 1.565a +2.194 



r = Exp 



s = 101.22 X Exp 



27r( ;/-(3-2a) 



(-7.301 - 1.565a) 




Defining u = BDql, the deflection angle a{6) can be rewritten as 



1 I ^ ^ps 

a = —am ' 



+ b + 0{u - Ups) 



that works out to 



a{a) 



-77 + 2.123- 1.565a +10.478 



1 



3 - 2a 



(87) 

(88) 
(89) 
(90) 

(91) 
(92) 



The values of the observables are tabulated below. We see that the values of a,b continue to 
increase from the Schwarzschild values (a = 1, 6 = —0.4009) as we increase the tidal charge. We also 
observe that the deflection angle a{a) increases from the Schwarzschild value as the tidal charge is 
increased as opposed to the decrease caused by ordinary scalar fields (e.g., JNW scalar field) [8]. 
This difference due to the tidal charge a is particularly manifest in the WH region corresponding to 
1 < a < |. For a ^ |, the deflection angle a{a) increases more than three times compared to the 
value a(|) for the Schwarzschild BH. Such behavior could be interpreted as a signature for a WH as 
well as effect of the extra dimension or tidal charge. The behavior of the observables r and s too are 
very different from the Schwazschild BH (or the JNW scalar field configuration) for different values 
of a, especially at a ~ |. 



Table I 
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Sch: 
7 = 1 


7 = 1.2 


7 = 1.5 


7 = 1.8 


7 = 2 


a = 0.8 


a = 1 


a = 1.2 


a = 1.4 


a 

(rad.) 


6.36 


6.53 


6.72 


6.87 


6.96 


6.58 


7.89 


10.63 


20.22 


(/i arcsec) 


16.87 


20.56 


26.02 


31.36 


35.00 


16.87 


16.87 


16.87 


16.87 


S 

{/i arcsec) 


0.0211 


0.0205 


0.0197 


0.0189 


0.0185 


0.0261 


0.0726 


0.3047 


3.1618 


(magn.) 


6.82 


6.82 


6.82 


6.82 


6.82 


6.59 


5.57 


4.31 


2.50 




2.59 


3.16 


4.00 


4.82 


5.38 


2.59 


2.59 


2.59 


2.59 


a 


1 


1 


1 


1 


1 


1.03 


1.22 


1.58 


2.74 


b 


-0.4009 


-0.4292 


-0.4692 


-0.5073 


-0.5321 


-0.4163 


-0.3895 


-0.0641 


1.6963 



V. Some features of the WH solutions 
A. Massive EMS WH 



The EMS solutions (38,39) correspond to an equation of state po + Pp + 2p_L = for the WH 
case 7 > 1 since po = Pp and —pp =pe = Pip- The equation is the hmiting case of the dark equation 
of state p = wp where w < — i. (The phantom, equation of state is more stringent as it requires 
w < —1 which is certainly not the case here.) The first observation is that the total asymptotic 
gravitating mass M = mj is positive. It can be calculated in various ways: by the ADM calculation 
[24] or from the Einstein energy complex or even directly from the Eddington-Robertson expansion 
of the centrally symmetric metric in isotropic coordinates [55] 



ds^ = ( 1 ^ + + 0{M^/p^) ] dt^ - 



2 



1 + hlM. + 5^ + O(MVp^) ) [dp' + p\de' + sin^ Od^')] (93) 
p Ip J 

Thus, the scalar field eff'ect is already contained in the metric functions A{p), B{p) in terms of 

M = Twy. This mass M is the gravitating mass and the test particles respond to it per se: there is 
in fact no way of measuring the bare m if a scalar field gravitationally couples to it. The Eddington- 
Robertson parameters for the Buchdahl solution (38) are ai = /3i = 71 = 1, and the post-PPN 
parameter (5i = | — The Buchdahl PPN parameters ai, 71 arc exactly the same as those in 
the Schwarzschild solution and at this level EMS theory is indistinguishable from it. However, the 
deviation appears at the post-PPN level and only finer and second order deflection measurements 
can reveal the value of Si. It is known that 5i = 1 (or, 7 = 1) corresponds to Schwarzschild solution 
while 61 ^ 1 would indicate a genuine deviation from it. The second order effect in deflection {albeit 
still in the weak fleld) can be easily calculated by using the metric (93) involving M and 5i and 
Eq.(6). 

However, due to the nonlinearity of the field equations, the total amount of WEC violating scalar 
matter flAWEC in spacetime is slightly different from — g as a result of the generalized Gauss theorem 
in curved spacetime. The exact difference can be seen from the volume AWEC integral, which, for 
two sides of the WH becomes 

^AWEC = 2x^^^y j j pn\/^sm6dpd6d(fi 

(94) 



- (\/f?T) " — (1 - ^) 
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As 7 1 (Schwarzschild case), f^AWEC ~* 0, as expected. Using the WH range, 7 > 1, wc sec that 
the Qawec < 0. To express ^Iawec in terms of scalar charge </, recall the two cases (i) and (ii) 
discussed in Sec. IIIA: If wc take k = —2, $ real or q ~ ra\J 7^ — 1, then ^awec — — 27^) ' 

q > {) while for n = +2, $ imaginary or g = im^J 7^ — 1, one has VIawec — ^ 27^) ■ 

integrate from pps to 00, we get VIawec — — and similarly the imaginary version. The 

positive correction term proportional to slightly diminishes the quantity ^awec away from the 
value — q, but this is due to nonlinear effects. For either of the two values of k and q, we see that 
= + q^ = m^7^. The situation is the following: the WH is attractive and test particles, 
after being pulled into the throat from both the mouths, attain zero acceleration there. They can 
re-emerge into the other mouth by maintaining extra outward radial accelerations from being pulled 
in again [56]. However, light rays are captured at p = Pps but rays that pass close to it suffer higher 
deflection angles than those due to Schwarzschild BH. 

B. Massless EMS WH 

The zero mass WH with its metric given by Eqs.(51) or (52) is a stable configuration (see Ref.[41] 
for details). The scalar field satisfies, to first order in $(~ ^), a sourceless equation -§^{p^^) = 0, 
yet the observers at a finite asymptotic p measures a flux 47rg' of the scalar charge though, in reality, 
there is no source. That is why the configuration is called a "charge without charge" . 

C. GMCFM I solution 



The GMCFM I solution is well discussed in the literature. Its asymptotic physical mass can be 
found from the Einstein complex as follows [56] 

M = LimP^^^ (95) 
p^o° V2AB 

which works out to M = m(x + A) which is finite and positive. The WH is sustained entirely by 
the negative pressures as po = 0. With the identifications A = 1 + e and >c = — e, we immediately 
find that M = m. Therefore, energetically, it is still like the same Schwarzschild spacetime while, 
kinematically, the null geodesies reveal that the strong field behavior is different from that of the 
Schwarzschild, as the observables obtained by Whisker [37] show. 

D. GMCFM II solution 



The GMCFM II solution exhibits certain remarkable features. Let us suspend the unit 2m = 1 
and restore m for better comparison. First, wc see that the constituent matter is that of dark 
radiation given by pD — {Px + 2p±) = 0. Second, it is impossible to ascertain the bulk effect directly 
from A{x) as it does not contain a. Therefore, we proceed as follows. The integration of the Einstein 
complex of energy gives the asymptotic physical mass of the solution [57] : 

M = i(2a + m) (96) 

which contains the bulk effect. From the Newtonian limit of gooix) = A(x) = 1 — — , the Kcplcrian 
mass is always m but the asymptotic mass M, except in the special case a = is different. This 
feature is unlike the Buchdahl solution where both masses are the same M (= rwy). When we 
interpret Eq.(96) as a relation of the type of Eq.(35) 

= (to^ + 4a^ + 4am) (97) 
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we see that, unlike the EMS case, there is an extra interaction term Aam between the mass and 
the Weyl charge contributing to the total mass M. We see that the Weyl charge is additive to m, 
unlike the scalar charge. Therefore, a has the dimension of mass. Expressing a in units of M such 
that a = eM, and using again the PPN expansion in standard coordinates for a central mass M: 



\ X ) 

(l + ^ + 1^ + ...) dx- - x\de^ + sin^ ^V) (98) 

and identifying the metric (63) with it, we find that a\= (i\ = 2(2 — e), 71 = 1,5\ = |(2£^ — 7£ + 8). 
At e = |, one recovers the Schwarzschild values ai = /?i = 71 = (5i = 1, as expected. When 
a = |M, we have M = m from Eq.(96). It turns out that first order tests can provide a value of 
e = I and can still be indistinguishable from GR, but (weak field) second order deflection would 
check if (5i = 1. If the measured value of 5i differs from unity, then the brane model of stars would 
stand as a possible contender to that in the EMS theory. 

The total amount of energy in spacetime due of bulk stress is 



^Totai = 2 X j j j PDy/-9 sin 9dxd9dip = m- ^ (99) 

which is negative for the WH range ^ < a < 3m. (Recall that Xth = a hes below Xpg = 3m). The 
integral is zero either for a = |m or, for a = |M since, in both cases, m = M. The AWEC by 
definition is 

/ \ /'2-7r PIT /.oo 

^AWEC = '^^ ( — I / / / pD\/—gsin9dxd9dip = m (100) 

\OTtJ Jo Jo Jxth=a 

which is independent of a\ This shows that the entire negative energy —^J is concentrated below 

the throat < a:; < a. However, there are some restrictions on the values of a. The pointwise WEC 
and the AWEC are satisfied here for a > For a < WEC is violated as pu < but Q.awec 
= —ioo which is unphysical. The energetics of the brane model requires further study but we see 
that the situation is very unlike the massive EMS WHs in which both WEC and AWEC arc violated 
for 7 > 1, the value of i^AWEC being proportional to the scalar charge —q, as shown before. 

Third, for a < we find that the gravitating mass M is decreased from the Schwarzschild 
value TO, the latter occurring at a = For a > we find that M > m, which suggests that 
the presence of the positive Weyl charge a strengthens the attractive force beyond that due to the 
Schwarzschild BH. This explains why there is an enhancement in the deflection angle. The two 
surfaces Xpg = 3to, and xth = ci coincide when a = 3to,. At this extremal situation, there occurs 
photon capture as the divergence in the deflection angle a{a) at a = 3to (which is the same as 
a = |) show in Eq.(92). At nonextremal situations that we have considered, the throat lies below 
Xps, and the light rays do not reach the throat. We are interested, as mentioned before, in the range 
^ < a < 3to and not in the range a > 3m as, in this case, it is the throat that covers the photon 
sphere, not the other way around. 

Finally, fourth, the WEC {po > 0) is locally preserved for a > ^ and so is AWEC which we 
saw to be independent of a. What about the ANEC violation? Let us consider the volume ANEC 
integral 

^ANEC = 2 X y J J {pD +Px)\/^sm0dxd0difi (101) 

2(.T-2TO)(2a-3TO) 



2am 



x'^{2x — 3to)2 



-s/^dx (102) 



16 



which contains only the radial pressure Px- (We have not included the transverse components of 
pressure as inequalities associated with p± refer only to ordinary matter [54]). That is, the WH 
is maintained by negative radial pressure [58] . Now, pointwise NEC violation po + Px < occurs 
when ^ < a < 3m, 2m < x < oo. Unfortunately, the ANEC integral diverges logarithmically on 
[a, oo). Such a divergence seems to be a generic feature of -R = traversable WHs [27]. However, it 
turns out that ft an EC < for a < x < N where N is any arbitrarily large but finite number. But 
as TV increases to infinity, so does ^Ianec though not as rapidly. Thus, to have a reasonable WH 
with finite amount of ANEC violating matter, one perforce needs to join the WH to the exterior 
vacuum Schwarzschild spacetime at a certain value of the coordinate radius x, as is actually done 
with self-dual B = WHs [54]. This peculiarity distinguishes the ANEC violation due to the bulk 
effect from the violation due to sign reversed scalar field. In the latter case, the ^^anec tapers off 
smoothly at the asymptotic region. 

VI. Summary 

Gravitational lensing by WHs in the strong field limit is a new possibility that has not been 
explored so far though pioneering theoretical works on lensing by black holes or naked singularities 
exist in the recent literature. On the other hand, WH solutions are physically important, many of 
their properties have been widely discussed and applied to interpret several outstanding problems in 
astrophysics [59-62]. It is thus only natural that their analytic lensing properties be investigated as 
well. Several static, spherically symmetric WH solutions are known, both in the EMS theory and in 
the brane theory. Some of the brane world solutions represent a synthesis of BH and WH spacetimes 
thereby providing a more advanced and richer premise for the strong field lensing analysis. We have 
undertaken a moderately c;oniprehensive investigation here. Certain intrinsic features of the lensing 
objects in question are also analyzed. 

The sign reversed kinetic term in Eq.(30) yields regular, symmetric WH solutions for the range of 
values 7 > 1 (The two options for k have been spelled out in the text.) For this range, the presence 
of the scalar charge increases the Schwarzschild mass m, that is, the ADM mass M ~ m'j > m. 
The WH throat surface lies inside the photon sphere which has been stipulated to play the limit of 
the strong field. That is, we have been considering situations in which light rays pass very close to 
the exterior of the photon sphere but obviously do not reach the throat. It was shown that massive 
WHs in the EMS theory (M 7^ 0) produce significantly different values of deflection angles, and 
other observables as tabulated in Table I. In contrast, in the case of naked singularity, $ real, k > 0, 
7 < 1, there is a decrease from the Schwarzschild mass, that is, M <m. This explains why, in this 
case, the deflection angles always show lesser values than those in the Schwarzschild case [8]. The 
strong field lensing results thus show that the EMS scalar field exerts stronger gravitational pull to 
light than that by the Schwarzschild BH. It should however be remarked that an increase in ADM 
mass does not generally imply an increase in the deflection angle. This follows from the fact that 
the bending features are unaffected by a conformal factor whereas a conformal factor does change 
the ADM mass. 

The massless WH (M = 7717 = 0) corresponds to 7 = but m ^ 0. This is just a mathematically 
admissible possibility. The spacetime is asymptotically flat at the two mouths. Lensing by these 
objects is interesting due to the fact that it can reveal the presence of the geometric curvature caused 
by the scalar field alone. It turns out that such configurations also possess a photon sphere and 
behave like ordinary deflectors. 

The GMCFM I solution is treated as a brane world BH in the literature [37], but it is actually 
a traversable WH. One recalls [52] that such class of solutions can be a BH only when it is trivially 
Schwarzschild. Otherwise, it is either a naked singularity or a WH. However, the throat radius Xth = 
2m is hidden under the photon sphere justifying the application of the strong field analysis. The 
GMCFM II solution has been investigated in detail here. Though the minimum impact parameter 
Ups is exactly the same as that in Schwarzschild BH, the spacetime itself is intrinsically very different. 
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Table I affords a comparison of the values of strong field observables for WHs with those of 
Schwarzschild BH. It shows that the separation (s) between the first image and the rest increases 
from that due to Schwarzschild BH with increasing 7 (> 1) but the increase is more spectacular 
in the brane world WHs, especially in the region a ^ 1.4. This suggests that the outermost image 
would be better visible in this case. The flux ratio (rm) or relative magnification always remains the 
same in the EMS WHs but is considerably more than those in the brane world WHs. These features 
are peculiar enough to observationally distinguish the lensing sources under consideration. 

Lensing phenomena in the WH environment offers a good possibility that one might detect the 
presence not only of a WH, which is by itself interesting, but also of the presence of naturally 
occurring exotic matter much advocated on galactic or cosmological scales. VLBI observations of a 
clean system that is devoid of intervention by accretion phenomena, can help us pick up the right 
model or at least set limits on the observables. Still, it is not unlikely that observations will favor 
just the usual Schwarzschild BH more than any other advanced generalized solution considered here. 
Again, we note that several astrophysical phenomena (like 7-ray bursts) can also be explained by 
invoking new inputs (like negative energy fields or exotic matter) [62], Thus, assuming that the 
center of our galaxy hosts, instead of a BH, a WH threaded by exotic matter, situated at a distance 
DoL = 8.5 kpc from the center of the Sun, then the angular position of the set of relativistic images 
in the limit n —>■ 00 would be 6*00 17 iJi arcsec. We have used this value as a basis for calculating 
s in Table I. However, due to considerable demagnification of relativistic images, one would need a 
resolution of the order of 0.01 fi arcsec and if this refinement is technologically attained in future, 
then the observational limits can either accommodate or rule out WH candidates as possible lensing 
sources. 

Acknowledgments 

We are indebted to Valerio Bozza for suggesting several improvements on an earlier version of 
the manuscript. We are thankful to Guzel N. Kutdusova and Guo Cheng for useful assistance. This 
work is financially supported in part by the TWAS-UNESCO program of ICTP, Italy, and by the 
Chinese Academy of Sciences and in part by the National Basic Research Program of China under 
Grant No. 2003CB716300 as weU as by the NNSFC under Grant No. 90403032. 

References 

[I] P. Schneider, J. Ehlers, and E.E. Falco, Gravitational Lenses (Springer- Verlag, Berlin, 1992). 
[2] S. FritteUi, T.P. Kling, and E.T. Newman, Phys. Rev. D 61, 064021 (2000). 

[3] K.S. Virbhadra and G.F.R. Ellis, Phys. Rev. D 62, 084003 (2000). 
[4] C.-M. Claudel, K.S. Virbhadra, and G.F.R. Ellis, J. Math. Phys. 42, 818 (2001). 
[5] V. Bozza, S. Capozziello, G. lovane, and G. Scarpetta, Gen. Relat. Grav. 33, 1535 (2001). 
[6] E.F. Eiroa, G.E. Romero, and D.F. Torres, Phys. Rev. D 66, 024010 (2002); For strong field 
retrolensing, see: E.F. Eiroa and D.F. Torres, Phys. Rev. D 69, 063004 (2004) 
[7] K.S. Virbhadra and G.F.R. Ellis, Phys. Rev. D 65, 103004 (2002). 
[8] V. Bozza, Phys. Rev. D 66, 103001 (2002). 
[9] A. Bhadra, Phys. Rev. D 67, 103009 (2003). 

[10] V. Bozza, F. De Luca, G. Scarpetta, and M. Sereno, Phys. Rev. D 72, 083003 (2005). 

[II] S.E. Vasquez and E.P. Esteban, Nuovo Cim. 119B, 489 (2004). 

[12] M.S. Morris and K.S. Thorne, Am. J. Phys. 56, 395 (1988); M.S. Morris, K.S. Thorne, and 
U. Yurtsever, Phys. Rev. Lett. 61, 1446 (1988). 

[13] S.M. CarroU, Living Rev. Rel. 4, 1 (2001); R.R. CaldweU, Phys. Lett. B 545, 23 (2002). 

[14] C. Charmousis, V. Onemh, Z. Qiu, and P. Sikivie, Phys. Rev. D 67, 103502 (2003); Z.-H. 
Zhu, Int. J. Mod. Phys. D 9, 591 (2000); C.-M. Claudel, Proc. Roy. Soc. Loud. A 456, 1455 
(2000). 

[15] V. Onemh, astro-ph/0510414 

[16] F.S.N. Lobo, Phys. Rev. D 71, 084011 (2005); S. Sushkov, Phys. Rev. D 71, 043520 (2005). 



18 



[17] J.G. Cramer, R.L. Forward, M.S. Morris, M. Visser, G. Benford, and G.A. Landis, Phys. 
Rev. D 51, 3117 (1995). 

[18] M. Safonova, D.F. Torres, and G.E. Romero, Phys. Rev. D 65, 023001 (2002). 
[19] S.J.M. Tejeiro and R.E.A. Larranaga, gr-qc/0505054 . 

[20] K.S. Virbhadra, D. Narasimha, and S.M. Chitre, Astron. Astrophys. 337, 1 (1998). 

[21] A.G. Agnese and M. La Camera, Phys. Rev. D 51, 2011 (1995). 

[22] K.K. Nandi, A. Islam, and James Evans, Phys. Rev. D 55, 2497 (1997). 

[23] L.A. Anchordoqui, S.P. Berghaffa, and D.F. Torres, Phys. Rev. D 55, 5226 (1997) 

[24] K.K. Nandi, B. Bhattacharjee, S.M.K. Alam, and J. Evans, Phys. Rev. D 57, 823 (1998). 

[25] A. Bhadra and K. K. Nandi, Mod. Phys. Lett. A 16, 2079 (2001). 

[26] A. Bhadra and K.K. Nandi, Phys. Rev. D 64, 087501 (2001). 

[27] K.K. Nandi, Yuan-Zhong Zhang, and K.B.Vijaya Kumar, Phys. Rev. D 70,127503 (2004). 
[28] K.K. Nandi, Yuan-Zhong Zhang, and K.B.Vijaya Kumar, Phys. Rev. D 70, 064018 (2004). 
[29] K.K. Nandi, Yuan-Zhong Zhang, Phys. Rev. D 70, 044040 (2004). 

[30] Arunava Bhadra, Ion Siniaciu, K.K. Nandi, and Yuan-Zhong Zhang, Phys. Rev. D 71, 
128501 (2005). 

[31] N. Dadhich, R. Maartens, P. Papadopoulos, and V. Rezania, Phys. Lett. B 487, 1 (2000). 

[32] R. Rochaand C.H. Coimbra-Araiijo, JCAP 0512, 009 (2005) (astro-ph/0510318, astro-ph/05093^); 
C.H. Coimbra-Araujo, R. Rocha, and I.L. Pedron, Int. J. Mod. Phys. D 14, 1883 (2005) 
(astro-ph/0505132). 

[33] L. RandaU and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999); ibid, 83, 4690 (1999). 

[34] R. Guedens, D. Clancy, and A.R. Liddle, Phys. Rev. D 66, 043513 (2002). 

[35] E.F. Eiroa, Phys. Rev. D 71, 083010 (2005). 

[36] A.S. Majumder and N. Mukherjee, astro-ph/0403405 . 

[37] R. Whisker, Phys. Rev. D 71, 064004 (2005). ' 

[38] D.F. Torres, E.F. Eiroa, and G.E. Romero, Mod. Phys. Lett. A 16, 1849 (2001); A 16, 973 
(2001). 

[39] A. Janis, E. Newman, and J. Winnicour, Phys. Rev. Lett. 20, 878 (1968); 

[40] M. Wyman, Phys. Rev. D 24, 839 (1985); R. Gautreau, Nuovo Cimento B 62, 360 (1969). 
For the equivalence of JNW, Wyman and Buchdahl solutions, see: A. Bhadra and K.K. Nandi, Int. 
J. Mod. Phys. A 16, 4543 (2001). 

[41] C. Armendariz-Pi'con, Phys. Rev. D 65, 104010 (2002). 

[42] H.A. Buchdahl, Phys. Rev. 115, 1325 (1959). 

[43] J.A. Wheeler, Phys. Rev. 48, 73 (1957). 

[44] H.G. Ellis, J. Math. Phys. 14, 104 (1973); 15, 520(E) (1974). For scattering problems 
in such spacetimes, see: G. Clement, Int. J. Theor. Phys. 23, 335 (1984); L. Chetouani and G. 
Clement, Ge. Relat. Gravit. 16, 111 (1984). 

[45] V. Perlick, Phys. Rev. D 69, 064017 (2004). 

[46] T. Shiromizu, K.I. Maeda, and M. Sasaki, Phys. Rev. D 62, 024012 (2000). 

[47] M. Visser and D.L. Wiltshire, Phys. Rev. D 67, 104004 (2003). 

[48] S.W. Hawking, G.T. Horowitz, and S.F. Ross, Phys. Rev. D 51, 4302 (1995). 

[49] C. Germani and R. Maartens, Phys. Rev. D 64, 124010 (2001). 

[50] R. Casadio, A. Fabbri, and L. Mazzacurati, Phys. Rev. D 65, 084040 (2002). 

[51] R. Casadio and B. Harms, Phys. Lett. B 487, 209 (2000). 

[52] K.A. Bronnikov, V.N. Melnikov, and H. Dehnen, Phys. Rev. D 68, 024025 (2002). 

[53] N. Dadhich, S. Kar, S. Mukherjee, and M. Visser, Phys. Rev. D 65, 064004 (2002). 

[54] M. Visser, S. Kar, and N. Dadhich, Phys. Rev. Lett. 90, 201102 (2003). 

[55] S. Weinberg, Gravitation and Cosmology, (John Wiley, New York, 1972), p. 183 

[56] L.H. Ford and T.A. Roman, Phys. Rev. D 53, 5496 (1996). 

[57] K.S. Virbhadra, Phys. Rev. D 60, 104041 (1999). 

[58] K.A. Bronnikov and S.-W. Kim , Phys. Rev. D 67, 064027 (2003). 



19 



[59] Ke-Jian Jin, Yuan-Zhong Zhang, and Zong-Hong Zhu, Phys. Lett. A 264, 335 (2000). 

[60] D. Hochberg and T.W. Kephart, Phys. Rev. Lett. 70, 2665 (1993); Phys. Lett. B 268, 377 
(1991); Gen. Relat. Gravit. 26, 219 (1994). 

[61] D.H. Coule, Phys. Lett. B 450, 48 (1999); D.H. Coule and K. Maeda, Class. Quant. Gravit. 
7, 955 (1990). 

[62] V. Prolov and I. Novikov, Phys. Rev. D 48, 1607 (1993). 

[63] F. Hoyle, G. Burbidgc, and J.V. NarHkar, Proc. Roy. Soc. Lond. A 448, 191 (1995); Mon. 
Not. Roy. Astron. Soc. 267, 1007 (1994). 



20 



